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Abstract 

We have investigated hadron screening masses, the chiral condensate, and the pion 
decay constant close to the deconfinement phase transition in the confined phase of QCD. 
The simulations were done in the quenched approximation, on a lattice of size 32"^ X 8. 
We examined temperatures ranging from 0.75Tc up to 0.92rc. We see no sign of a temper- 
ature dependence in the chiral condensate or the meson properties, but some temperature 
dependence for the nucleon screening mass is not excluded. 



1 Introduction 



Despite extensive studies of the QCD phase transition [1-3] and the phases on either 
side of it, httle is known with certainty about the changes in the excitation spectrum. 
In particular, the properties of the (quasi)-particle spectrum in the confined phase close 
to Tc require a more thorough understanding, as the temperature dependence of hadron 
masses and other hadronic parameters will lead to observable consequences in current, 
and especially forthcoming, heavy ion collisions. For example, a temperature dependence 
of the p mass and width leads to large modifications of the dilepton spectrum [4,5]. 

The temperature dependence of various hadronic properties has been addressed using 
a number of different approaches. These calculations yield different predictions for the be- 
haviour of some quantities. For example, in chiral perturbation theory it has been found 
that the chiral condensate and decrease, while and mjv increase [6] with increas- 
ing temperature. In an approach which uses a Skyrme-type Lagrangian, a relationship 
between various hadronic quantities was obtained [7] which suggests that all quantities 
should behave in the same way with increasing temperature. On the other hand, QCD 
sum rule predictions (eg. [5,8,9]) yield a decreasing nucleon mass and only a slow variation 
of nip. On the other hand, it has been argued [10] that hadron masses calculated from 
finite temperature sum rules are temperature independent up to 0{T^). 

Since the above approaches have limited applicability at high temperatures, one may 
hope that a lattice calculation of the variation with temperature of any of these quantities 
below the phase transition may shed some light on these discrepancies. 

At finite temperature the notion of the mass is somewhat ambiguous. One has to 
distinguish the pole mass defined through the pole position of the propagator from the 
screening mass obtained from the hadronic correlation functions at large spatial separa- 
tions [11-13]. These two will coincide only if the ordinary zero temperature dispersion 
relations also hold at non-zero temperature. Most of the lattice calculations at finite 
temperature which have determined screening masses [12-20], concentrated on the region 
T > Tc. The results are compatible with the propagation of deconfined quarks. The few 
calculations performed for T < were, however, all at rather small temperatures, and 
no significant changes from the T = masses were seen [21]. 

In this paper we present results from a study of hadronic properties in quenched QCD 
at temperatures between 0.75rc and 0.92rc. We will discuss the temperature dependence 
of the pion, rho and nucleon screening masses, the chiral condensate, and the pion decay 
constant. Details of the simulations performed here are presented in section 2. Section 3 
contains the results pertaining to the chiral sector, and section 4 those for the nucleon 
and /9-meson. Finally, in section 5 we present our conclusions. 

2 Simulations 

Previous studies have shown that the quenched approximation yields, for both the zero 
temperature masses and the high temperature screening masses, results that are within 
10% of those obtained using full QCD with both two and four fermion fiavours. Fur- 
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thermore, one sees a rapid, simultaneous change in the Polyakov loop and the chiral 
condensate in both the quenched approximation and in full QCD. This suggests that the 
hadronic properties are largely dependent on the gluon sector. We have thus chosen to 
use the quenched approximation for this study. 

We have to point out however, that close to the critical temperature the fluctuations 
associated with a second order phase transition may have a large effect on the temperature 
dependence of the screening masses. In this region the two flavour case may thus yield 
different results to both the four flavour and the quenched cases. The studies of the 
chiral transition in two flavour QCD [22] suggest, however, that these effects only become 
relevant very close to Tc. Hence we expect that our present analysis is not yet affected by 
these effects. 

We used a lattice of size 32"^ X 8 at coupling constants of (3 = 5.90, 5.95, 5.975 and 
6.0, giving a range in temperature from 0.75rc to 0.92rc. We have set the temperature 
scale using the calculation of the non-perturbative beta function by the TARO collab- 
oration [23], which indicates that the critical coupling lies at /3 ~ 6.06. This value is 
somewhat larger than those quoted in [24-26]. It is, however, supported by a direct 
determination of the critical point on 32"^ X 8 lattices, which is currently in progress [27]. 

Studies of the temperature dependence of hadronic properties close to Tc need large 
spatial volumes in order to prevent flips between coexisting phases. The volume used here 
was seen to be large enough to prevent flips to the deconfined phase, even at (3 = 6.00. 

We used four staggered valence quarks with masses niq = 0.05, 0.025 and 0.01 for the 
point sources, and a mass of niq = 0.01 for the wall sources. It is known that quenching 
introduces chiral logarithms which may spoil the chiral limit. However, they do not seem 
to present a difficulty for the quantities and within the parameter range studied here [28]. 

Details of the configurations generated are shown in Table 1. The first set of con- 
figurations was generated using a Kennedy-Pendleton heat bath algorithm (indicated by 
HB), and the second set using a combination of heat bath and over-relaxed algorithms 
(indicated by OR). For the pure heat bath algorithm we used 50 sweeps per configuration, 
whilst for the combined algorithm we used four over-relaxed sweeps per heat bath, and 
40 heat bath sweeps per configuration. 



/5 


T/Tc 


No. Confs 


niq 


5.90 


0.75 


150 (HB) 


0.05, 0.035, 0.025, 0.01 


5.95 


0.85 


150 (HB) 


0.05, 0.025, 0.01 


5.975 


0.90 


150 (HB) + 100 (OR) 


0.05, 0.025, 0.01 


6.0 


0.92 


350 (OR) 


0.05, 0.025, 0.01 



Table 1: The number of quenched configurations generated at each coupling is listed 
above, along with the valence quark masses used. The algorithm used was either pure 
heat-bath (HB) or a combined heat-bath and over-relaxation algorithm (OR). 

The correlators obtained for each channel have been fitted using both the full corre- 
lation matrix as well as only the diagonal elements. Errors have been calculated using 
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both the jackknife procedure and the standard + 1 definition. We found that in most 
cases both fitting procedures yield almost identical values for the fitted parameters and 
the jackknife errors, in particular when the errors on the propagators were small or when 
the fit interval could be limited to large distances from the source. Based on this, and the 
findings of [29], we have quoted masses based on fits using the diagonal elements of the 
correlation matrix for those cases where the fit with the full matrix did not work. The 
errors refiect the error from jackknife estimates, plus an estimate of the additional error 
coming from varying the fitting range. 



Within the context of QCD sum rules as well as chiral perturbation theory the variation 
of hadronic parameters is closely related to the temperature dependence of fermionic and 
gluonic condensates. We thus will first discuss an analysis of the temperature dependence 
of the chiral condensate. It is known to change drastically at the deconfinement transition 
point in quenched QCD. However, a detailed investigation of its temperature dependence 
below Tc on large lattices did not exist so far. 

We have calculated the chiral condensate directly from the trace of the fermion matrix 
((?/'V')se) and from the pseudoscalar (PS) and scalar (SC) propagators, which corrects for 
the linear dependence on the quark mass [30]. The latter method is given by the formula 



The meson propagators on the right hand side are taken at zero four momentum. 

Figure 1(a) shows the value of the chiral condensate as a function of bare quark mass 
at a coupling (3 = 6.0, corresponding to 0.92rc. One can see quite clearly in the figure 
that both methods extrapolate to the same value at zero quark mass. The same holds 
true at the other couplings studied. 

The temperature dependence of the chiral condensate is shown in figure 1(b), where we 
plot the ratio of the finite temperature chiral condensate extrapolated to zero quark mass 
with that at r = 0. For comparison results for two [22], three [31-33] and four [34,35] 
fiavour QCD, and another quenched study [21] have been included in the figure. These 
additional values have not been extrapolated to zero quark mass, which may lead to 
an underestimate of the temperature effects. However, it is clear from the figure that 
the chiral condensate does not change significantly until one is actually at the critical 
temperature. Fven at the largest beta value, corresponding to T = 0.92rc, the drop in 
the condensate is not statistically significant. 

A second quantity of considerable interest for understanding the temperature depen- 
dence of the chiral sector of QCD is the pion decay constant. It can be determined directly 
from the relevant matrix element [30] 



3 Chiral Sector 




(3.1) 
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Figure 1: Figure (a) shows the chiral condensate as a function of the bare quark mass, 
at (3 = 6.00. The upper curve is obtained using the value taken directly from the trace 
of the fermion matrix, the lower curve after subtracting the linear dependence on the 
quark mass. The T = results have been taken from [37]. Figure (b) shows the finite 
temperature chiral condensate normalised to its value at zero temperature. The circles 
represent our values, the inverted triangles a quenched calculation in [21] at Nf = 4, both 
extrapolated to vanishing quark mass and normalised to the data of [37,38]. The triangles 
show the result using the intercept method for a four fiavour calculation [34,35] at Nf = 8, 
squares a 2 + 1 fiavour calculation at niq = 0.0125 in [31-33] at Nf = 4 and stars the 
results for 2 fiavours at niq = 0.02. [22]. 
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Figure 2: Figure (a) shows the square of the pion screening mass at /3 = 6.00 plotted as 
a function of the bare quark mass. The circles represent our data, the triangles the data 
of [36] at zero temperature. Figure (b) shows at T > normalised to the value at 
r = [37], at a quark mass of = 0.01 in lattice units. For reference we include the 
values obtained for the matrix element above the phase transition as well. 



At zero temperature is related to the chiral condensate and the pion mass through the 
Gell-Mann, Oakes, Renner (GMOR) relation, 

^Ifn = ?7^q(V'V')(mq=0) • (3.3) 

This relation is expected to be valid as long as chiral symmetry remains spontaneously 
broken. 

A calculation of and a test of the GMOR relation at finite temperature requires a 
calculation of the pion mass. We have calculated at four values of the temperature the 
quark mass dependence of the pion screening mass. In all cases we find the expected 
behaviour for a Goldstone particle, = A^m^. The amplitude, A^r, does not show any 
temperature dependence. In figure 2(a) we have plotted the pion screening mass squared 
as a function of quark mass at a coupling of (3 = 6.0, along with the results at zero 
temperature. It is clear that both sets of points lie on the same straight line. There is no 
dependence on temperature; the relation oc niq holds true, up to T = 0.92rc. 

The temperature independence of the amplitude, A^, together with our finding of 
a temperature independent chiral condensate up to T = 0.92rc indicates that the pion 
decay constant will also be temperature independent up to this temperature, if the GMOR 
relates still holds. Indeed we have found, at all couplings, that the pion decay constant 
determined from eqn. 3.2 and extrapolated to zero quark mass agrees very well with that 
extracted from the GMOR relation using the value of the chiral condensate extrapolated to 
zero quark mass. Thus the GMOR is indeed valid at temperatures as high as T = 0.92rc. 

In figure 2(b) we have plotted the value of the pion decay constant at the lightest 
quark mass used, niq = 0.01, normalised to the value at zero temperature as a function 
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of T/Tc. As expected we see no temperature dependence up to T = 0.92rc 



4 Nucleon and meson masses 



The calculation of liadron masses usually proceeds through more or less involved fits 
to the hadron correlation functions Cyi{z) and requires a discussion of the fitting range 
etc. Another more direct way to discuss the long distance behaviour of the correlation 
functions is the analysis of local screening masses obtained from ratios of the hadron 
correlation functions at neighbouring lattice sites. 

These are obtained by solving the following equation for the mass m: 

Cnjz) ^ cosh(m(z - NJ2)) 
Cn{z + 1) cosh{m{z + 1 - NJ2))' ^ ' ^ 

For large distances these ratios stabilize and yield the lowest energy in the quantum 
number channel selected by the correlation function under consideration. As we are 
considering correlation functions at large spatial separation, the relation of these energies 
to screening masses differs at finite temperature for bosonic and fermionic states. The 
latter also receive a contribution from the non-vanishing Matsubara energy, po = ttT, i.e. 

E^ = ml + k sin^ {tt/K) . (4.2) 
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0.02 0.04 0.06 0.08 

Figure 3: The rho screening mass at /3 = 6.00 plotted as a function of the bare quark mass. 
The circles represent our results using point sources, except for niq = 0.01 where the value 
obtained from wall sources is given. The triangles represent the zero temperature results 
from [36,37,39,40]. 

The rho screening mass at T = 0.92rc is shown in figure 3 as a function of the quark 
mass, including the zero temperature values. The local screening masses as a function of 
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distance, as well as lines indicating the value obtained from a fit to the full propagator 
here and at zero temperature are shown in figure 4(a). It is clear from both figures that 
the rho screening mass does not show any significant dependence on temperature. As 
can be read off Table 2, where we have summarized our results for the screening masses, 
the rho screening mass is, within errors, identical in both pseudo-vector and vector-tensor 
channels, indicating that the continuum flavour symmetry is restored in the vector meson 
spectrum. The table also includes our numbers for the scalar, the /o, where as usual 
the quarkline disconnected contribution has been omitted. Note that independent of the 
temperature, the /o exhibits a constant ratio to the p mass, nifg/nip ~ 1.2(1). 

Finally, we turn to the nucleon. This was examined using wall sources at a quark mass 
of 0.01 in lattice units. The local screening masses for T = 0.92rc, and the screening mass 
obtained from a full flt, are shown in flgure 4(b) along with the result obtained at zero 
temperature. The increase in temperature clearly has a dramatic effect on the nucleon. 
However, most of this can be understood purely in terms of the fact that the lowest 
momentum of the nucleon is not zero, but ttT, as discussed above. 

The screening mass at T = 0.92rc, extracted using eqn. 4.2 with the assumption that 
A; = 1, indicates that the mass rises slightly: mN(0.92rc) = (1.1 it 0.03)mN(r = 0). We 
note, however, that there is also the possibility of a modiflcation of the energy dispersion 
relation at flnite temperature, which may lead to deviation of k from unity [41]. The error 
given above is from the statistical error alone, as we cannot determine the systematic error 
which may be introduced by our assuming that eq. 4.2 with A; = 1 is applicable. 
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Figure 4: Figure 4(a) shows the local rho screening mass at T = 0.92rc, and niq = 0.01 
in lattice units, from wall sources. The solid and dotted lines show the fltted value and 
errors, the long dashed lines the zero temperature values from [39]. Figure 4(b) shows 
the local nucleon screening mass at T = 0.92rc, and niq = 0.01 in lattice units, from wall 
sources. The solid and dotted lines show the fltted value and errors, the long dashed lines 
the zero temperature values from 
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r article 


™ n ni (\M\ 

filq — U.Ui 1^ VV J 


filq — U.UZOI^r J 


?Tiq — U.UtjOl^r j 


™ n DKii'p^ 

?Tiq — U.UOI^r j 






a k; on 






TT 


0.267(2) 


0.411(3) 


0.484(2) 


0.572(2) 


Pl 


0.61(2) 


0.71(3) 


0.75(2) 


0.80(2) 


P2 


0.55(4) 


0.63(3) 


0.70(3) 




ai 


0.85(10) 


1.0(1) 


1.0(1) 




fo 


0.67(5) 


0.77(3) 


0.83(2) 




TT 


0.35(5) 


0.52(6) 


0.60(4) 




N 


0.95(5) 








(3 = 5.95 


TT 


0.255(5) 


0.391(3) 




0.550(3) 


pl 


0.50(3) 


0.58(2) 




0.72(2) 


P2 


0.54(3) 


0.63(3) 




0.70(3) 


ai 


0.8(1) 


0.8(1) 




1.0(1) 


fo 


0.62(3) 


0.71(2) 




0.83(2) 


TT 


0.34(3) 


0.46(4) 




0.62(3) 


N 


0.95(5) 








(3 = 5.975 


TT 


0.255(2) 


0.397(2) 




0.552(2) 


pl 


0.58(6) 


0.63(3) 




0.72(1) 


P2 




0.65(3) 




0.70(2) 


ai 








0.85(10) 


fo 


0.63(6) 


0.74(3) 




0.84(1) 


TT 


0.29(3) 


0.44(2) 




0.60(1) 


(3 = 6.00 


TT 


0.247(2) 


0.383(2) 




0.538(2) 


pl 


0.49(1) 


0.57(2) 




0.70(1) 


P2 


0.47(1) 


0.55(2) 




0.68(1) 


ai 


0.72(3) 


0.85(10) 




0.97(7) 


fo 


0.58(2) 


0.71(3) 




0.83(1) 


TT 


0.30(2) 


0.41(2) 




0.59(1) 


N 


0.86(3) 









Table 2: Masses obtained from fits to the full propagator. W denotes wall sources, P 
denotes point sources. Details of the fitting procedure are discussed in section 2 
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5 Conclusions 



We have analyzed the temperature dependence of various hadronic parameters. Probably 
the most unexpected result is that the most basic observable, the chiral condensate, does 
not show any significant temperature dependence up to T = 0.92rc. In view of this our 
observation that there is no sign of temperature dependence in the pion decay constant, 
pion screening mass, /o screening mass or the rho screening mass may not come as a 
surprise. There is some, albeit inconclusive, indication of a temperature dependence for 
the screening mass of the nucleon. 

Our calculations have been performed in the quenched approximation of QCD. The 
details of the temperature dependence may change in the case of QCD with two fiavours, 
where the transition is expected to be second order rather than first order as in the pure 
SU (3) gauge theory. However, as discussed in Section 2 we do not expect these differences 
to show up for temperatures less than O.QTc. 
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